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A comparison of Clover and Wilson spectroscopy in the presence of 
dynamical quarks 

S. Collins, R.G. Edwards, U.M. Heller, and J. Sloan'* 
^SCRI, Florida State University 

We present preliminary results of light hadron spectroscopy using valence, tadpole-improved. Clover fermions 
on an ensemble of gauge configurations generated with 2 flavors of staggered fermions at /3 = 5.6. We compare 
the slope and intercept of the curve Mv vs. Mpg for Clover and Wilson fermions. We show that a higher order 
chiral perturbation theory ansatz works very well for chiral extrapolations. 



1. INTRODUCTION 

In this talk, we present preliminary results from 
a detailed comparison of Clover and Wilson light 
spectroscopy at a moderate inverse lattice spacing 
(about 2 GeV). We discuss two analytic methods; 
1) (tadpole-improved) Symanzik improvement of 
the fermion action and 2) higher order chiral ex- 
trapolations. In addition, we use correlated fits 
to multiple smearing functions and an automated 
plateau finder with different selection criteria in 
an attempt to reduce systematic fitting errors. 



2. THE SIMULATION 

The gauge configurations were generated by 
the HEMCGC collaboration [ with /? = 5.6 
and 2 flavors of staggered fermions at a mass of 
am = .01. Inversions were performed on 100 con- 
figurations separated by 20 trajectories, fixed to 
Coulomb gauge. We use tadpole improvement in 
our calculations; uq = .867 is obtained from the 
plaquette. 

Quark propagators were inverted for each k us- 
ing a gaussian source; at some k's a wall source 
was also used. Mixed k hadrons were formed 
between propagators generated in the same run. 
For the V and PS mesons, each spatial (mesonic) 
smearing yields two meson interpolating fields; 
one with spin structure F and the other with spin 
structure F * 79. Similarly, we measure two dif- 
ferent A+ states; Sz = 3/2 and Sz — 1/2. 



3. FITTING 

A technique we find very useful is to perform a 
correlated multi-state fit to multiple correlation 
functions which have the same quantum num- 
bers. Two types of fit are to a "vector" of 
smeared source and local sink and to a "matrix" 
of smeared source and sink correlations. Even if 
only a single energy is used in the fit ansatz, mul- 
tiple smearings can help eliminate the problem 
of spurious plateaus because the multiple prop- 
agators will have different overlaps with excited 
states. A notation we find useful is to call a single 
state fit involving two smearings a "lc2s fit." 

In choosing our fits, we have tried to avoid 
human intervention by using automated plateau 
finders. I discuss two of our criteria here; "C", 
which tends to pick very short plateaus (con- 
servative) and " D" , which prefers long plateaus 
(aggressive). Because correlation functions are 
noisier at later times, the trade-off between dif- 
ferent criteria is between statistical uncertainties 
and systematic fitting errors; an aggressive cri- 
terion will give smaller error bars but may have 
more excited state contamination in the ground 
state mass. The goal is to use the most ag- 
gressive criterion for which the statistical un- 
certainty dominates fitting errors. The defini- 
tion of C is mm{T„iin + /^prop)', for D it is 
max((5 * Ndof /SM), where Nprop is the number 
of propagators being fit, Q is the confidence level 
of the fit, and SM is the fit uncertainty in the 
mass. 

To perform chiral fits, we generated 201 boot- 
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strap ensembles of gauge configurations, with cor- 
responding sets of propagators. After choosing a 
plateau for a particular fit, we repeated the fit on 
each bootstrap ensemble, resulting in an ensemble 
of 20 i mass and amplitude values. Note that, un- 
like [ 1^ , the plateau was kept fixed over all boot- 
strap ensembles. We would like to implement the 
method of [ |^ , which allows the plateau to vary 
between bootstrap ensembles, as this will move 
systematic fitting errors into bootstrap fluctua- 
tions where they can be statistically quantified. 

Note that the mass ensembles depend both on 
the type of fit (i.e. Icis vs. Ic2s) and on the fit 
criterion used (i.e. C vs. D). A too aggressive cri- 
terion results in large in the chiral fit; this is 
because excited state effects (systematic errors) 
are about the same size as the statistical errors. 
This allows investigation of systematic fitting er- 
rors; our results should be consistent between dif- 
ferent criteria and types of fits. 

4. Kc 

Our first example of a higher order chiral ex- 
trapolation is in the determination of Kc- In 
Fig. 1^, we show Mpg vs. 1/k for the clover ac- 
tion, along with quadratic and linear fits. The 
quadratic fit includes all six points shown, while a 
linear fit is possible only for the first three points. 
The values of Kc extracted from these fits are 
.i40809(i6) and .i40782(i6), respectively A typ- 
ical problem when determining Kc is that it "runs 
away" ; linear extrapolations to Kc change as 
lighter k's are included. This does not seem to be 
the case for the quadratic fit; omitting the lightest 
two k's from the fit results in a Kc of .i4079(2), 
which is consistent with the full quadratic fit. Be- 
cause of this we trust the quadratic fit; even our 
lightest K values are probably too heavy to be in 
the linear regime. 

5. MESONS 

In Fig. H we present our results for My vs. 
Mpg for the Wilson and clover actions; both 
degenerate and non-degenerate n combinations 
are included. The vertical dotted lines repre- 
sent the p, K* and (j) masses. The p and K* 
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Figure I. Linear and quadratic fits of Mpg vs. 
1/kc for clover action. 

were chosen by requiring that the ratio My / Mps 
match experiment. The pseudoscalar mass corre- 
sponding to the 4> was chosen using the formula 
Af^j = 2m\r, — m^. The mass values for both 
V and PS were obtained from lc2s vector fits us- 
ing criterion D; the two smearings were the same 
spatially and differed by 70 in spin structure. 
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Figure 2. Quadratic and cubic chiral fits to My 
for clover and Wilson actions. 

Chiral fits were performed using four different 
ansatze; 

2- 1 : Co + CiMps^ 

3- 1 : Co + CaMps' + CsMps^ 
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4-2 : Co + C2Mps'^ +CiMps^ 
4-1 : Co + C2Mps^ + CsMps^ + CaMps^ 

where the first number indicates the order of the 
polynomial and the second is the number of co- 
efficients set to zero. The 2-1 (quadratic) ansatz 
is just lowest order xPT; the 4-2 (quartic) ansatz 
introduces a higher order term that goes like . 
The cubic term in the 3-1 (cubic) and 4-1 (full 
quartic) ansatze were motivated by the one loop 
xPT correction to baryon masses; at this confer- 
ence we learned of [ ||], in which one loop Mps'^ 
corrections to the vector meson mass are derived. 

The results of chiral fits were qualitatively sim- 
ilar for Wilson and clover. In both cases, a 
quadratic fit breaks down after the sixth point 
is included, while both the cubic and quartic (4- 
2) fits were able to include many more points 
(typically 11-18 points). The cubic fit, however, 
seemed to do a better job of describing data at 
even higher masses which could not be fit; it and 
the quadratic fits are shown in Fig. ||. The full 
quartic (4-1) ansatz only allowed one or two more 
points to be fit than either a cubic or quartic 
ansatz; this is a sign that there is no signal to 
resolve the additional parameter. The numer- 
ical results of our fits are presented in Tab. ^ 
We have also performed fits using only degener- 
ate K combinations and obtain consistent results. 
Note that we can resolve the difference between 
Mp at Mt^ = .ISA/p and at M,r = 0; we obtain 
aMp = .341(4), .400(4) for Wilson and clover, re- 
spectively. 



Table 1 

Chiral fit results for My vs. Mps- indicates 
light k's were omitted. 



Fit 


Wilson 


Clover 


Co 




Co 




2-1 


.337(6) 




.395(6) 




3-1 


.336(4) 


1.23(2) 


.394(4) 


1.10(3) 


4-2 


.337(6) 


1.13(5) 


.403(6) 


0.93(2) 


3-1* 


.334(4) 


1.25(2) 


.393(4) 


1.10(3) 



Note that the cubic fits are consistent with and 
have smaller error bars than the quadratic fits. 



This indicates that the higher mass points, which 
have smaller fiuctuations, are controlling the fit, 
i.e. the lightest (and most expensive) k's are too 
noisy to contribute useful information. To test 
this, we omitted all k combinations lighter than 
the third degenerate k (6 points for Wilson and 4 
for clover) and repeated the 3-1 and 2-1 fits. The 
quadratic ansatz only fit one degree of freedom 
(3 points) for clover and none for Wilson. The 
cubic fits, however, succeeded to the same Umax 
and gave results with the same error bars. This 
indicates that it was unnecessary to simulate be- 
low the strange quark mass in order to calculate 
the physical p mass. 

6. BARYONS 

For baryons, we have not yet analyzed our Wil- 
son results. In addition, we have results with wall 
smearings at only two of our six lightest k's. We 
will correct these deficiencies in the near future. 




Figure 3. Comparison between fit criteria and 
between single and multi-propagator fits for the 
K — .1375 nucleon. Vertical lines indicate Tmin 
chosen by criteria D and C (D is always leftmost). 
Upper right plot is superposition of effective mass 
for shell and wall sources; upper left plot is mass 
values for Ids fit as function of Tmm- 

One interesting result for the nucleon was that 
we were unable to obtain good chiral fits when 
criterion D was used for Ids propagator fits. We 
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interpret this as an example of too aggressive a 
fit criterion leading to statistically significant ex- 
cited state contamination of the ground state en- 
ergy. This can be corrected by either using a less 
aggressive fit criterion (i.e. C), or by including a 
second smearing (i.e. wall). This can be seen in 
Fig. H the Ids criterion D fit uses T„iin = 1! This 
is changed to 5 and 7 for the C criterion and lc2s 
fit, respectively. In Fig. ^ we show that using cri- 
terion C results in successful chiral fits, but that 
C is probably too conservative. The A behaves 
similarly; criterion D gives mass values which can 
not be fit, while criterion C works well. As with 
the vector, a cubic ansatz allows much heavier 
quark masses to be included. 
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Figure 4. Chiral fits to Mnucieon generated from 
fcls fits and criterion C. 



7. CONCLUSIONS 

Our most striking result is that a cubic ansatz 
for chiral extrapolations of the vector mesons 
works so well, and over such a large mass range. 
Our results seem to indicate that the two k values 
lighter than Ustrange were unnecessary, since their 
omission from fits caused negligible changes. This 
is especially important given that the light k's are 
both the most expensive and the most sensitive 
to finite volume effects. A potential drawback to 
this approach was our inability to clearly distin- 
guish between cubic and quartic ansatze, and the 



lack of a strong signal to justify including both 
terms at the same time. 

Table 2 

Comparison of dimensionless ratios 

Observable Wilson Clover Expt. 

Td .425(6) .444(5) l99~ 

{Mn- M^)/Mk' - .39(5) .491 

Mn/MK' .575(11) .536 



Our comparison between Wilson and clover ac- 
tions has only been completed for mesons; the 
results were reasonable. For both the intercept 
and slope of My vs. Mpg, we saw 10-20% dis- 
crepancy. This is consistent with a naive esti- 
mate of 500 MeV * a « .25 for Wilson discretiza- 
tion errors, and suggests that the O(a^) clover 
discretization errors are roughly 5% at this lat- 
tice spacing. A comparison between clover re- 
sults and experiment is less encouraging; as can 
be seen from Tab. |[ the ratio Jd = Mk* * {Mk' — 
Mp)/{^K - ^i) is 10% too small and the ratio 
(Mn - Ma) /Mk- is 20% too small. We hope that 
this effect is due to either finite volume or incor- 
rect dynamical fermion content; if this is not the 
case then either clover is not correctly removing 
0(a) errors or there is some source of error which 
we do not understand. 
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